Antiferromagnetic order in multi-band Hubbard models for iron-pnictides 
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We investigate multi-band Hubbard models for the three iron 3d-t2g bands and the two iron 3d-eg 
bands in LaOFeAs by means of the Gutzwiller variational theory. Our analysis of the paramagnetic 
ground state shows that neither Hartree-Fock mean-field theories nor effective spin models describe 
these systems adequately. In contrast to Hartree-Fock-type approaches, the Gutzwiller theory 
predicts that antiferromagnetic order requires substantial values of the local Hund's-rule exchange 
interaction. For the three-band model, the antiferromagnetic moment fits experimental data for a 
broad range of interaction parameters. However, for the more appropriate five-band model, the iron 
Eg electrons polarize the t2g electrons and they substantially contribute to the ordered moment. 
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Since their recent discovery, the iron-based high-Tc su- 
perconductors have attracted tremendous attention both 
in theory and experiment. From a theoretical point of 
view, these systems are of particular interest because 
their conduction electrons are less correlated than those 
of other high-Tc superconductors. In contrast to the 
cuprates, the pnictides' undoped parent compounds are 
antiferromagnetic metals at low temperatures, not insu- 
lators. However, the electronic mass is enhanced by a 
factor of two which indicates that electronic correlations 
are quite substantial in the pnictides, too. 

The theoretical description of the pnictides' normal 
phase turned out to be a difficult problem. Standard 
density-functional theory (DFT) grossly overestimates 
the size of their magnetic moment in the antiferromag- 
netic ground state. For example, in LaOFeAs experiment 
finds a staggered moment of m = (0.4 . . . 0.8)/iB [IHl] 
whereas DFT calculations predict moments of m w 1.8/xb 
or larger 0, [1]. For other pnictide compounds, the com- 
parison is equally unfavorable. 

Perturbative many-electron theories for the pnictides 
face the problem that all five iron Sd-bands contribute 
significantly to the band structure near the Fermi energy. 
Even when the effective Hubbard interaction U between 
pairs of electrons on the same iron atom is smaller than 
the effective 3d bandwidth W, U < W/2, atomic charge 
fluctuations with ionization degree AA^ > U/W are im- 
probable in the ground state. The description of itinerant 
antiferromagnetism in such correlated multi-band sys- 
tems is beyond perturbative mean-field approximations 
such as Hartree-Fock (HF) theory. In studies based on 
dynamical mean- field theory (DMFT), only the param- 
agnetic phases of these models could be investigated so 
far, see e.g., Ref. @. 

In this work, we employ the Gutzwiller variational the- 
ory (GT) which describes the ground state and quasi- 
particle excitations of multi-band Fermi liquids; for an 
application to nickel, see Refs. 0, As compared to 
DMFT, our method is numerically much less costly. This 



enables us to resolve the small energy differences between 
the paramagnetic and antiferromagnetic phases in the 
pnictides (Neel temperature Tn ~ 140 K). From our 
study of a five-band model for the iron 3d-t2g and 3c?- 
Cg bands we find that the paramagnetic phase is much 
more stable than predicted by DFT or HF calculations. 
Only for large values of the local exchange interaction, 
an antiferromagnetic ground state is found; in a related 
study [Si], only paramagnetic states were examined. 

We investigate the two-dimensional five-band Hubbard 
model 

^= E <'l'4MS-,fc'..+E^c..^^o+^ioc, (1) 
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where Hq describes the transfer of electrons with spin 
a =t,i between iron atoms on lattice sites i,j in the 
3d-{eg,t2g) orbitals 5,6'. The transfer parameters t'^'j 
are taken from Ref. [l3|- The bare bandwidth of the 
3d electrons is W = 4.8 eV, and there are, on average, 
six electrons on every iron atom. The local Hamiltonian 
Hc,i describes the Coulomb interaction of the iron elec- 
trons. Frequently the following form for the Hubbard 
interaction is employed, 

b;cr b(^)b' 
a. a' 

H^^ = J{b, b') (4,t^i^,c,, ^c,,^ + h.c.) (2) 
fc(#)fc' 

+ Jib,b')clJ,,,c,,c,.^. 
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Here, we dropped the lattice-site indices and introduced 
the abbreviations f =i, I =tj and Ua,a'{b, b') = U{b, b')- 
Sa-,c7' J{b,b'), where U{b,b') and J{b,b') are the local 
Coulomb and exchange interactions. 
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Even in cubic symmetry, however, the Hamiltonian ([2]) 



is incomplete. The full Hamiltonian reads Hq = iJ, 



(1) 



E^^^ with 



(2) 



^ 5(t,t';i",e)c,Vc: 



h.c. 



Here, t and e are indices for the three t2g orbitals with 
symmetries xy, xz, and yz, and the two Cg orbitals with 
symmetries u = 3z^—r^ and v — x^ — y^. The parameters 
in ([3]) are of the same order of magnitude as the exchange 
interactions J(&, b') and, hence, there is no a-priori reason 
to neglect them. Of all the parameters U{b,b'), J{b,b'), 
A{t), T{t), S{t, t'; t", e) only ten are independent in cubic 
symmetry. In a spherical approximation, they are all 
determined, e.g., by the three Racah parameters A, B, C. 
For details on the multiplet structure of d-shells we refer 
to Sugano's textbook [llj. 

In this work, we work with the orbital averages J cx 
J2b^b' "^(^'^Oj ^^'^ (X X]f,5:^f)' U{b,b') of the exchange 
and the inter-orbital Coulomb interaction. They are re- 
lated to the intra-orbital interaction U = U{b, b) via 
U' = U ~ 2J . Due to this symmetry relation, the three 
values of U, U' , and J do not determine the Racah pa- 
rameters A, B, C uniquely. Therefore, we further assume 
that the approximate atomic relation C/B = 4 is satisfied 
in solids, too. In this way, the three Racah parameters 
and, consequently, all parameters in Hq are functions 
of U and J. This permits a meaningful comparison of 
our results with those of previous work where, as an ad- 
ditional approximation, the parameters U(b,b'), J(b,b') 
were assumed to be orbital independent, see, e.g., Ref. Q. 
We note that the Gutzwiller method is applicable to the 
general form ([1]) of the atomic Hamiltonian. 

We approximate the true ground state of -ff in ([T]) by 
the Gutzwiller variational wave function 



(4) 



where |^o} is a single-particle product state, and the local 
'Gutzwiller correlator' is defined as 



p. = 5]Ar|r},,(r| 



(5) 



Here, we introduced variational parameters Ar for each 
of the atomic multiplet states |r)i, i.e., the eigenstates 
of Hc,i- Note that the single-particle state Iv^o) is also 
a variational object which we determine from the mini- 
mization of the variational energy functional that results 



from the wave functions ([4]). The Hartree-Fock theory is 
a special case of the Gutzwiller theory which we obtain 
by setting Ar = 1 for all |r). 

The evaluation of expectation values with respect 
to dU poses a difficult many-particle problem. As 



shown in Refs. 12, it can be solved exactly in the 
limit of infinite spatial dimensions. The analytic en- 
ergy functional derived in this limit can be used as an 
approximation for finite dimensional systems. On top 
of the approximate ground-state description which is 
provided by the Gutzwiller wave function one can also 
calculate quasi-particle band structures for a compari- 
son with data from angle-resolved photoemission spec- 
troscopy (ARPES) Q. 

The Gutzwiller energy functional contains two differ- 
ent sets of variational parameters. The local multiplet 
occupations are governed by the parameters Ar in the 
correlation operator ([5]). As shown in Refs. [13|, the 
single-particle wave function |5'o) is the ground state of 
an effective single-particle Hamiltonian 
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and 



with renormalized electron transfer parameters ^ 

variational parameters rjl'^ which govern the orbital and 
spin dependent local densities. There are N\ = 2^" = 
1024 parameters Ar while the number of parameters r/^'^ 
is much smaller, 7V,,.p — 5 for the paramagnetic case and 
A^rj.af = 10 for the antiferromagnetic case. Nevertheless, 
the minimization of the energy functional with respect 

to the parameters ry^'^ is numerically expensive because 
any change of these parameters results in a minimization 
cycle with a full momentum-space integration. In order 
to obtain the required energy resolution, these integrals 
have been calculated on a momentum-space grid with up 
to 3 • 10*^ triangles in the Brillouin zone. 

In order to test the reliability of our approach for the 
five-band model, we first compare our results for the par- 
tial densities with those from paramagnetic DMFT cal- 
culations. In Fig. [1] we show the density of electrons in 
each orbital as a function of U for fixed ratio U/J = A. 
The full symbols give the GT result for the simplified lo- 
cal Hamiltonian ([2]), = Hq + H^^}''; open symbols give 
the DMFT results. Obviously, the agreement between 
the GT and DMFT is very good despite the fact that the 
particle number is not perfectly conserved in the DMFT 
calculations 
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Fig. [T] shows a common feature of multi-band model 
systems. The local Coulomb interaction induces a sub- 
stantial charge flow between the bands because, for the 
local Coulomb interaction, it is energetically more favor- 
able to distribute electrons equally among the bands. 
However, the bands described by Hq are usually ex- 
tracted from a DFT calculation whose predictions for 
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FIG. 1. (Color online) Orbital densities in GT (full symbols) 
and in DMFT (open symbols) as a function of U (with U / J = 
4) for the simplified local Hamiltonian He ~ Hq \ see ((2|. 



FIG. 3. (Color online) Magnetic moments in the five-band 
model in Hartree-Fock (HFT) and Gutzwiller theory (GT) 
as a function of J for [/ = 2.2 eV, 2.6 eV, 3.0 eV. 



the Fermi surface reproduce experimental data reason- 
ably well. Therefore, we argue that the artificial charge 
flow as seen in Fig.[T]is a consequence of the double count- 
ing of Coulomb interactions. Since the (paramagnetic) 
Fermi surface found in DFT reproduces its experimen- 
tally determined shape, we assume that the same holds 
for the paramagnetic orbital densities. For each value 
of the interaction parameters we therefore choose orbital 
on-site energies eb = which lead to a paramagnetic 
ground state with the same orbital densities as in DFT. 
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FIG. 2. (Color online) Local charge distribution (a) and spin 
distribution (b) for the paramagnetic optimized Gutzwiller 
wave function for U = 2.6 eV, J = 0.4 eV (left columns) and 
for the Hartree-Fock wave function (right columns). 

Next, we investigate the local charge distribution and 
spin distribution in the paramagnet. In Fig. [^^ we show 
the probabilities to find an atom with < n < 10 elec- 
trons in I^-qP') for U = 2.6 eV and J = 0.4 eV. As 
seen from Fig. [2^, local charge fluctuations with A7V = 
|iV — 6| > 2 are essentially forbidden. Hartree-Fock- 
type mean-field approximations cannot describe this ef- 
fect properly. The probabilities to find the local spins 
< s < 5/2 in l^-g") are shown in in Fig. Wp- The 
broad spin distribution in Gutzwiller theory is very sim- 
ilar in Hartree-Fock theory and shows that the system 
is far from the local-moment situation where we would 
solely find atoms with Hund's-rule spin s ~ 2. Con- 
sequently, the elementary excitations of the five-band 



Hamiltonian ([T} with parameters from Ref. [lo| and 
[U = 2.6 eV, J = 0.4 eV) are not well approximated by 
an effective spin model. 

We now turn to the magnetic properties of our five- 
band model III). As seen in Fig. [31 Hartree-Fock the- 
ory predicts an antiferromagnetic phase for moderate 
2.2 eV <U< 3.0 eV for all J > 0. The HF moments 
are strongly orbital dependent. For small values of J, 
these moment are anti-aligned which leads to a phase 
with a small total moment. At critical values J = Jc{U), 
the HF moments align which results in a large-moment 
phase. Similar phases have been reported in LDA+U 
calculations 16]. 

Our Gutzwiller theory shows that the correlated para- 
magnetic state is stable over a wide range of Coulomb 
and exchange-interaction parameters, e.g., J < 0.4 eV 
for 2.2 eV < U < 3.0 eV in Fig. H The small- J phase 
with orbital-dependent moments as seen in Hartree-Fock 
theory is absent. The robustness of the paramagnetic 
phase against symmetry breaking is in agreement with 
earlier findings [17] that Hartree-Fock theory generally 
overrates the importance of orbitally ordered phases in 
the ground state of multi-band systems. 

In Gutzwiller theory, the transition to an antiferromag- 
net with a large moment, m > 2/iB, is as abrupt as in 
Hartree-Fock theory. As seen in Fig. [31 we find an anti- 
ferromagnetic state with a small moment m = 0.4/iB only 
in a small region of parameter space, e.g., U = 2.2 eV and 
0.42 eV < J < 0.45 eV. Since this region in the U-J pa- 
rameter space is fairly small, we do not consider it very 
likely that the iron pnictides fall into this parameter re- 
gion. The explanation for their small antiferromagnetic 
moment should have a more natural explanation. 

We have extended our analysis to a broader class of 
variational wave functions where we included the mixing 
of atomic configurations in the Gutzwiller correlator ([5]) 
of the form Ar,r' |r)ii(r'| for F 7^ F'. In order to keep the 
problem numerically tractable, we considered the dom- 
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inant 500 couplings for which |(«'^P*|r)ii(r'|\E'°P*)|2 is 
largest in the paramagnetic Gutzwiller state I^'q^*}- The 
inclusion of these additional variational parameters low- 
ered the energies of the paramagnetic and the antifer- 
romagnetic optimal states by almost equal amounts so 
that the phase diagram does not change noticeably. In 
particular, the region in phase space with a low magnetic 
moment increases only marginally. 



^4 U=2.2eV 

• • U=2.6eV 

• « U=3.0eV 



30 0.35 0.40 0.45 0.50 0.55 0.60 0.55 0.70 0.75 

J/eV 



FIG. 4. (Color online) Magnetic moments for the three-band 
model as a function of J for U = 2.2 eV, 2.6 eV, 3.0 eV. Full 
symbols: GT with the full local interaction He = H^^ ((2)); 
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Our findings for the five-band model contrast those of a 
recent study of a three-band model for the pnictides [18 1 
where the authors report a small magnetic moment over 
a large region of Coulomb and exchange interactions. In 
that study, the two Brf-e^-bands and their interactions 
H^^ with the 3d-i2£rbands, eq. were not considered. 
Moreover, the spin-flip and pair-exchange terms in the 
local Hamiltonian, the terms in eq. have been 
neglected. 

In Fig. m we show the Gutzwiller prediction for the 
magnetic moments of the three-band model for 2.2 eV < 
U < 3.0 eV with the full local interaction, He — 
Hq , and, for comparison, the result with He 
Note that we used the electron transfer amplitudes of 
Ref. [3l for the parameterization of Hq in ([T]) , and that 
we use a slightly different minimization algorithm. In 
contrast to the full five-band model, a large region with 
a small magnetic moment is found for the three-band 
model. When the spin-flip and pair-transfer terms H^ 
are included in the local interaction, the transition in the 
three-band model occurs at larger values of J and, for 
given J > Jc, the magnetic moment is generally smaller 
for the full local interaction than for density interactions 
only. We obtain qualitatively similar results when we 
switch off the coupling between the t2g bands and the 
Eg bands in Hq of our five-band model with parameters 
for Hq taken from Ref. [IQ] ('3-|-2'-band model'). 

The comparison of the results for the three-band model 



and the five-band model shows that the Coulomb cou- 
pling between the 3(i-t2g-bands and the 3(i-eg-bands is 
very important. In the five-band model the 3d-e(,-bands 
contribute substantially to the magnetic moment and 
they strongly spin polarize the 3d-t2g bands. In the 
three-band model (or the '3-|-2-band model'), the Sd-Cg 
electrons are essentially paramagnetic. Since the pre- 
dictions of the three-band model fit experiment on the 
magnetic moment of the pnictides whereas the five-band 
model fails, it remains to justify the assumption that the 
3d-eg electrons are irrelevant for the magnetism of the 
pnictides. In principle, the effective Coulomb parameters 
which couple the i2£(-electrons and the Cg-electrons could 
be much smaller than those used in our study. In such a 
'3-f 2-band model', the Sd-Cg electrons would not carry a 
magnetic moment or polarize the 3d-t2g-electrons. More 
important for the pnictides is the hybridization of the 
irons' 3d-electrons with the arsenic 4p electrons. In fact, 
DFT calculations locate the arsenic 4p-bands not too far 
from the Fermi energy. Therefore, a full 3(i-4p model 
may provide a satisfactory description of the paramag- 
netic and antiferromagnetic phases of the iron pnictides. 
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